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We investigate random walks on complex networks and derive an exact expression for the mean 
first passage time (MFPT) between two nodes. We introduce for each node the random walk 
centrality C, which is the ratio between its coordination number and a characteristic relaxation 
time, and show that it determines essentially the MFPT. The centrality of a node determines the 
relative speed by which a node can receive and spread information over the network in a random 
process. Numerical simulations of an ensemble of random walkers moving on paradigmatic network 
models confirm this analytical prediction. 
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From biology over computer science to sociology the 
world is abundant in networks. Watts and Strogatz J^J 
demonstrated that many of these real-world networks 
pose characteristic features like the small world and clus- 
tering property. Since classical network models do not 
display these features, new models have been developed 
to understand the structure underlying real-world net- 
works 1^13. The important finding is that there exists 
a scale-free (SF) network characterized by a power-law 
degree distribution that is also a characteristic for 
the World- Wide- Web (WWW) 0, and for many other 
networks in various disciplines 

The SF network has a heterogeneous structure. For 
instance, the WWW analyzed in Ref. Q includes nodes 
with degrees ranging from 0{1) to 0(10'^). The het- 
erogeneity leads to intriguing properties of SF networks. 
In the context of percolation, SF networks are stable 
against random removal of nodes while they are fragile 
under intentional attacks targeting on nodes with high 
degree 0, @ . Statistical mechanical systems on SF net- 
works also display interesting phase transitions 0, H, ■ 
In a transport process, each node in SF networks does not 
contribute equally likely. The importance of each node in 
such a process is measured with the betweenness central- 
ity 01 , which has a broad power-law distribution ■ 

In this paper, we study a random walk on general net- 
works with a particular attention to SF networks. The 
random walk is a fundamental dynamic process |l2j |. It 
is theoretically interesting to study how the structural 
heterogeneity affects the nature of the diffusive and re- 
laxation dynamics of the random walk ■ Those issues 
will be studied further elsewhere 0|. The random walk 
is also interesting since it could be a mechanism of trans- 
port and search on networks 0,0,0]. Those processes 
would be optimal if one follows the shortest path be- 
tween two nodes under considerations. Among all paths 
connecting two nodes, the shortest path is given by the 
one with the smallest number of links 0| . However the 



shortest path can be found only after global connectivity 
is known at each node, which is improbable in practice. 
The random walk becomes important in the extreme op- 
posite case where only local connectivity is known at each 
node. We also suggest that the random walk is a useful 
tool in studying the structure of networks. 

In the context of transport and search, the mean first 
passage time (MFPT) is an important characteristic of 
the random walk. We will derive an exact formula for 
the MFPT of a random walker from one node i to an- 
other node j, which will be denoted by (Ty), in arbitrary 
networks. In the optimal process it is just given by the 
number of links in the shortest path between two nodes, 
and both motions to one direction and to the other di- 
rection are symmetric. However, a random walk motion 
from i to j is not symmetric with the motion in the oppo- 
site direction. The asymmetry is characterized with the 
difference in the MFPT's. It is revealed that the differ- 
ence is determined by a potential-like quantity which will 
be called the random walk centrality (RWC). The RWC 
links the structural heterogeneity to the asymmetry in 
dynamics. It also describes centralization of information 
wandering over networks. 

We consider an arbitrary finite network (or graph) 
which consists of nodes i = 1, . . . , iV and links connect- 
ing them. We assume that the network is connected (i.e. 
there is a path between each pair of nodes oth- 
erwise we simply consider each component separately. 
The connectivity is represented by the adjacency ma- 
trix A whose element Aij = 1 (0) if there is a link 
from i to j (we set An = conventionally). In the 
present work, we restrict ourselves to an undirected net- 
work, namely Aij — Aji. The degree, the number of 
connected neighbors, of a node i is denoted by Ki and 
given by Ki = Y,j Aij. 

The stochastic process in discrete time that we study 
is a random walk on this network described by a master 
equation. The transition probabilities are defined by the 
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following rule: A walker at node i and time t selects one 
of its Ki neighbors with equal probability to which it hops 
at time t + 1, thus the transition probability from node 
i to node j is Aij /Ki . Suppose the walker starts at 
node i at time t = 0, then the master equation for the 
probability Pij to find the walker at node j at time t is 
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The largest eigenvalue of the corresponding time evolu- 
tion operator is 1 corresponding to the stationary dis- 
tribution Pj^ ~ limt^Qc Pij{t), i.e. the infinite time 
limit 21]. An explicit expression for the transition prob- 
ability Pij {t) to go from node i to node j in t steps follows 
by iterating Eq. Q 
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Comparing the expressions for Pij and Pji one sees im- 
mediately that 
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This is a direct consequence of the undirectedness of the 
network. For the stationary solution, Eq. (PJ implies that 
KiP°° = KjP°°, and therefore one obtains 

K 

P°° = — (4) 

with TV = J^i^i- Note that the stationary distribution 
is, up to normalization, equal to the degree of the node 
i — the more links a node has to other nodes in the 
network, the more often it will be visited by a random 
walker. 

How fast is the random walk motion? To answer to this 
question, we study the MFPT. The first-passage proba- 
bility Fij (t) from i to j after t steps satisfies the relation 



Inserting this series into Eq. (jBJ and expanding it as a 
power series in s, we obtain that 
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for i ~ i 
for j ^ i 
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A similar expression is derived in Ref. J^ll for the MFPT 
of the random walk in periodic lattices. 

It is very interesting to note that the average return 
time {Tii) does not depend on the details of the global 
structure of the network. It is determined only by the 
total number of links and the degree of the node. Since 
it is inversely proportional to the degree, the heterogene- 
ity in connectivity is well reflected in this quantity. In 
a SF network with degree distribution P{K) ^ K""', 
the MFPT to the origin To also follows a power-law dis- 
tribution P{To) - To'^'^^''\ The MFPT to the origin 
distributes uniformly in the special case with 7 = 2. 

Random walk motions between two nodes are asym- 
metric. The difference between (T^) and {Tji) for i ^ j 
can be written as (using Eq. ©) 
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where the last term vanishes due to Eq. Therefore 
we obtain 



(Tv) - (Trd = C;' - C-\ 



where Q is defined as 
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where P°° = Ki/Af and the characteristic relaxation 
time Ti of the node i is given by 
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The Kronecker delta symbol insures the initial condi- 
tion Pij{0) = 6ij {Fij{0) is set to zero). Introducing 
the Laplace transform f{s) = X^t^o ^~**-/^(^)' © be- 
comes Pij{s) ~ Sij + Fij{s)Pjj{s), and one has 

F,^^s)^{P,{s)-S.,)/P,,{s) . (6) 

In finite networks the random walk is recurrent so 
the MFPT is given by (T,,) = EZQtF^li^) = --^^'^(O)- 

Since all moments R\f = J^Zo^" {P^3W " ^3°°} of 
the exponentially decaying relaxation part of Pij (t) are 

finite, one can expand Pij as a series in s as 
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We call Ci the random walk centrality since it quantifies 
how central a node i is located regarding its potential 
to receive informations randomly diffusing over the net- 
work. To be more precise: Consider two nodes i and j 
with Ci > Cj. Assume that each of them launches a 
signal simultaneously, which is wandering over the net- 
work. Based on Eq. lO, one expects that the node with 
larger RWC will receive the signal emitted by its partner 
earlier. Hence, the RWC can be regarded as a measure 
for effectiveness in communication between nodes. In a 
homogeneous network with translational symmetry, all 
nodes have the same value of the RWC. On the other 
hand, in a heterogeneous network the RWC has a dis- 
tribution, which leads to the asymmetry in the random 
dynamic process. 
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FIG. 1: (a) r vs. K and (b) C vs. K calculated in the 
Barabasi-Albert network with N = 10000 and m = 2. The 
straight line in (b) has the slope 1. 



The RWC is determined by the degree K and r. The 
order of magnitude of the characteristic relaxation time 
r is related to the second largest eigenvalue {nota bene 
[2lj |) of the time evolution operator in Q: Pii{t) = 

P°° + T,\=2 4^''b[^^^\^ where a^^) and b'-^^ are the left 
and right eigenvectors, respectively, of the time evolu- 
tion operator belonging to the eigenvalue A^. If we order 
the eigenvalues according to the modulus (IA2I > IA3I > 
. . . > IAatI) the asymptotic behavior is Pii{t) — P°° ~ 

af^6f^A*2 and n « af Vl In IA2II. Thus the relax- 
ation time Ti has a node dependence only through the 
weight factor, which is presumably weak. On the other 
hand, the degree dependence is explicit. 

We examined the distribution of the RWC in the 
Barabasi-Albert (BA) network j22j. This is a model for 
a growing SF network; at each time step, a new node is 
added creating m links with other nodes which are se- 
lected with the probability proportional to their degree. 
We grew the network, solved the master equation nu- 
merically with the initial condition Pi{t = 0) = Sik, and 
calculated the relaxation time for each k. Figuren](a) 
shows the plot of r vs. K in the BA network of = 10'' 
nodes grown with the parameter m = 2. The degree is 
distributed broadly over the range 2 < K 400. On 
the other hand, the relaxation time turns out to be dis- 
tributed very narrowly within the range 1 < r ^ 2. We 
also studied BA networks of different sizes, but did not 
find any significant broadening of the distribution of r. 
So the RWC distribution is mainly determined by the 
degree distribution. In Fig. ^ (b) we show the plot of C 
vs. K in the same BA network. It shows that the RWC 
is roughly proportional to the degree. Note however that 
the RWC is not increasing monotonically with the degree 
due to the fluctuation of t as seen in Fig. ^ (a). 

The RWC is useful when one compares the random 
walk motions between two nodes, e.g., i and j with Ci > 




FIG. 2: (Color online) Time evolution of the fraction of walk- 
ers n that pass through a node as a function of (from top 
to bottom) the node index i, the node degree K and the 
RWC C of the BA network (left column) and the hierarchi- 
cal network (right column). The value of n at each time t is 
represented in the gray scale/color code depicted at the right 
border of each plot. 



Cj . On average a random walker starting at j arrives at 
i before another walker starting at i arrives at j. Now 
consider an intermediate node fc, which may be visited by 
both random walkers. Since (Tij) > (Tji), it is likely that 
a random walker starting at node k will arrive at node 
i earlier than at node j. Although this argument is not 
exact since we neglected the time spent on the journey to 
the intermediate node, it indicates that nodes with larger 
RWC may be typically visited earlier than nodes with 
smaller RWC by the random walker. If we interpret the 
random walker as an information messenger, nodes with 
larger RWC are more efficient in receiving information 
than nodes with smaller RWC. 

We performed numerical simulations to study the re- 
lation between the RWC and this efficiency. To quantify 
it, we consider a situation where initially all nodes in a 
network are occupied by different random walkers. They 
start to move at time t — 0, and we measure n^, the 
fraction of walkers which have passed through the node 
i, as a function of time t. It is assumed that the walkers 
do not interact with each other. They may be regarded 
as a messenger delivering an information to each node it 
visits. Then, with the information distribution uniformly 
initially, rii is proportional to the amount of information 
acquired by each node. The argument in the previous 
paragraph suggests that typically nodes with larger val- 
ues of RWC have larger value of rii at any given time. 
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The BA network l22l a nd the hierarchical network of 
Ravasz and Barabasi |23j were considered in the simula- 
tions. The hierarchical network is a deterministic net- 
work growing via iteration; at each iteration the net- 
work is multiplied by a factor M. The emergent net- 
work is scale- free when M > 3. Since it is a deter- 
ministic network, several structural properties are known 
exactly ■ We measured Ui in the BA network with 
m = 2 and N = 512 nodes and in the hierarchical net- 
work with M = 5 and TV = S'' nodes for < i < 2048, 
which are presented in the left and the right column of 
Fig. 13 respectively. The value of Ui is color-coded ac- 
cording to the reference shown in Fig. El 

The time evolution of rii is presented in three differ- 
ent ways. In the first row, the nodes are arranged in 
ascending order of the node index i. In the BA network, 
the node index corresponds to the time step at which 
the node is added to the network. The indexing scheme 
for the hierarchical network is explained in Ref. [2^. In 
the second row, the nodes are arranged in descending 
order of the degree K and in the third row they are ar- 
ranged in descending order of the RWC C. At a given 
time t, the plot in the first row shows that n is non- 
monotonous and very irregular as a function of the node 
index. As a function of the degree it becomes smooth, 
but still non-monotonic tendencies remain. However, as 
a function of the RWC, it becomes much smoother and 
almost monotonous. We calculated for each node i the 
time Tj' at which rii becomes greater than 1/2. In the BA 
network, among all node pairs satisfying < rj, 

only 3% violate the relation Ci > Cj, whereas the num- 
ber of pairs that violate the relation Ki > Kj is five times 
larger. 

In summary, we studied the random walk processes in 
complex networks. We derive an exact expression for the 
mean first passage time (see Eq. (O). The MFPT's be- 
tween two nodes differ for the two directions in general 
heterogeneous networks. We have shown that this differ- 
ence is determined by the random walk centrality C de- 
fined in Eq. H1U|) . Among random walk motions between 
two nodes, the walk to the node with larger value of C 
is faster than the other. Furthermore, it is argued that 
in a given time interval nodes with larger values of C are 
visited by more random walkers which were distributed 
uniformly initially. We confirmed this by numerical sim- 
ulations on the BA and the hierarchical network. One 
may regard the random walkers as informations diffus- 
ing through the network. Our results imply that infor- 
mation does not distribute uniformly in heterogeneous 
networks; the information is centralized to nodes with 
larger values of C. The nodes with high values of C have 
the advantage of being aware of new information earlier 
than other nodes. On the other hand, it also implies 
that such nodes are heavily loaded within an informa- 
tion distribution or transport process. If the network 
has a finite capacity, the heavily loaded nodes may cause 



congestions [I3 • Therefore much care should be taken of 
the nodes with high C values in network management. In 
the current work, we consider the random walks on undi- 
rected networks. The generalization to directed networks 
would be interesting. And in order to study congestion, 
the random walk motions with many interacting random 
walkers would also be interesting. We leave such gener- 
alizations to a future work. 
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